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ABSTRACT 

We present VLA Hi imaging data for a field in the NW of the galaxy cluster Abell 1367 
(z = 0.02) in an attempt to probe the effect environment has on the interstellar medium of 
_G late-type spiral galaxies. Several galaxies, like CGCG 097-087, show pronounced tails and 

Qh asymmetries, and 7 out of 10 show significant, several kpc offsets between the HI centroid and 

q the optical. We compare our results against a sample of optically bright, late-type galaxies 

5_i (spirals) across the central 1 .5 Mpc of the cluster taken from the Arecibo Galaxy Environment 

Survey (AGES). We find that these late-type spirals are predominantly found in the northern 
CS half of the cluster, especially those that are relatively gas-rich. We calculate the H I deficiency 

and find that the expected global trend for the H I deficiency of these spirals to increase with 
! projected proximity to the cluster core, seen in clusters like Coma and Virgo, is not observed. 

We classified the spirals into four evolutionary states, with the galaxies in each state sharing a 
similar degree of H I deficiency and optical colour. The common characteristics of the spirals 
C^) in each evolutionary state suggests they have been subject to similar environmental processes. 

Q\ Many of the spirals in the most common evolutionary state (moderate H I deficiency and blue 

C*~) colour) have an H I intensity maximum which is displaced relative to its optical counterpart. 

The orientation of these offsets and magnitude of their H I deficiencies together with data 
from other wavelengths provide observational evidence in support of varying degrees of ram 
c^~) pressure stripping and tidal interaction. In general, our results indicate that the H I disks of 

y—( bright late-type galaxies in the central part of the cluster are subject to both gas loss and 

morphological disturbance as a result of their interaction with the cluster environment. This 
. £^ provides further observational evidence of a more complex environment in Abell 1367 as 

compared to Virgo and Coma. 

Key words: galaxies: evolution — galaxies: ISM — galaxies: clusters: individual: 
(Abell 1367) 



1 INTRODUCTION 

For several decades it has been known that the fraction of spi- 
ral galaxies decreases when moving from a field population to a 
galaxy cluster core, with a corresponding increase in the fraction 
of early type (E+S0) galaxies; this is known as the morphology- 
density relation (Dressier 1980; Oemler 1974). Whether this re- 
lation originates during galaxy formation or is mainly the result 
of environment remains one of the fundamental questions in ob- 
servational cosmology. While previous work (Boselli & Gavazzi 
2006a; |Dressler et al.|1999[ |Poggianti et al.|1999| |Dressler|2004| 



* E-mail: t.c.scott@herts.ac.uk (TS) 



and references therein) suggests that the environment associated 
with clusters only plays a secondary role in producing lenticu- 
lar galaxies from spirals, there is ample evidence for the effect of 
cluster environment on galaxy evolution, such as the observed in- 
crease in the fraction of H I-deficient spiral galaxies toward clus- 
ter cores (e.g., Boselli & Gavazzi 2006a; [van Gorkom|2004l and 
references therein), where the Hi deficiency is defined as the log 
of the difference between the expected and observed Hi mass 
(see Table [5] table footnote h). There is strong evidence that in- 
falling cluster galaxies encountering a hot X-ray emitting intraclus- 
ter medium (ICM) for the first time experience rapid evolution as 
a result (Gavazzi et al. 2001b; [Sun & Murray|2002| |Kenney et al.| 
|2004| |van G orkom 2004 1. In these cases the interaction removes 
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the interstellar medium (ISM) from the galaxies (e.g., |Gavazzi| 
|et al.|1995||Bravo-Alfaro et al.|2000||Soranes et al.|2001| >. For in- 
falling spirals the gas motions resulting from the interaction are 
expected to lead to quenching of star formation. Although modest 
and short timescale (10 s years) increases in star formation are pre- 
dicted in specific circumstances (Fujita & Nagashima 1999, Kro- 
nberger et aL 2008 ). However the predicted ram pressure induced 
star formation was not observed in a study by Iglesias-Paramo et al. 
420041. 



Different types of mechanisms are proposed to explain how 
the cluster environment impacts evolution of its constituent galax- 
ies. The first type, which is most evident in spirals, consists of hy- 
drodynamic mechanisms arising from gas-phase interactions be- 
tween the cluster's hot ICM and a galaxy's ISM; these include ram 
pressure stripping ( |Gunn & Gott|1972|, viscous stripping l |Nulsen| 
|1982[>, thermal evaporat ion (Cowie & Songaila 1977) and starva- 
tion jBekki et al.|2002| |Fujita & Goto|2004| >. The second type is 
caused by gravitational effects produced either by close encounters 
with neighbouring galaxies (merging or tidal stripping) or by re- 
peated high speed galaxy-galaxy encounters (harassment; Moore 
|etal.|1996|[l999) . Also in this cate gory belongs gravitational inter- 
action between a galaxy and the cluster potential as a whole (Bekki 
[T99^[Natarajan et al.|2002) >. 

Hi in spirals provides an excellent tool to investigate these 
mechanisms as it is sensitive to both gravitational and hydrody- 
namic interactions (Dickey & Gavazzi 1991). There is conflicting 
evidence for and against a correlation between the fraction of H I 
deficient galaxies and the cluster X-ray luminosity, Lx, (Boselli 
& Gavazzi 2006a; |Giov"anelli & Haynes|1985[|Solanes et al.|200T) . 
However modelling by |Tonnesen| j2007) favours ICM stripping 
over gravitational effects as the primary cause of Hi deficiency in 
nearby cluster spirals. Modelling of a generic spiral infalling to a 
cluster predicts the removal of almost all of its H I during the first 
transit of the cluster core as a result of hydrodynamic stripping. 
During the transit a combination of high ICM density and maxi- 
mum orbital velocity hugely increases ram pressure stripping effi- 
ciency (Vollmer et al 2001a; Roediger & Briiggen 2007). The few 
individually modelled spirals in Virgo, e.g., NGC 4522 ( |Kenney| 



etaTT2 004; Voll mer et al.|2 008), NGC 4569 (Boselli et al.2006b; 
Vollmer et al. 2004} and in Coma, NGC4848 (Vollmer et al. 2001b), 



confirm that the ram pressure mechanism is operating. There are 
cases in Virgo, e.g., NGC 4654 (|Vollmer|2003[), NGC 4438 l|Boselli| 
|et al.|2005) , NGC 4254 ( |Chyzy|2008^ and (|Chung et al.|2007) as 
well as examples in other clusters ( Sun & Vikhlinin 2005 ; Moran 
|et al.|20 07 1 indicating that tidal mechanisms are operating in addi- 
tion to hydrodynamic mechanisms. 

The issue of galaxy evolution near clusters has become richer 
and more complex in recent years, since modelling suggests that 
under the hierarchical large-scale structure formation scenario, 
clusters are expected to grow predominantly by accretion of groups 
of galaxies rather than individual galaxies(Blummenthal et al. 
1984; Springel et al. 2001, although for an alternative model see 
Berrier et al. 2009) . Infalling groups, having lower velocity disper- 
sions than the clusters, provide opportunities for tidal interactions 
which may represent an important preprocessing stage in cluster 
galaxy evolution l|Mihos|2004||Fujita & Goto|2004||Dressler|2004|l. 



Wilman et aTT]l |2009| ) report strong evidence that transformation of 
spirals to SOs in clusters occurs preferentially in such groups as a 
result of tidal interactions at z ~0.4. Their results are based on anal- 
ysis of HST-ACS images of 179 spectroscopically confirmed group 
and 111 field galaxies from the Canadian Network for Cosmology 
(CNOC2) sample. This raises the question of what processes have 



perpetuated the transformation of late-type galaxies to SOs since 
then ( |Fasano et al.|2000") , given in the current epoch hydrodynamic 
mechanisms appear to be the predominant interaction type, but the 
most likely mechanisms for conversion of spirals to SOs are widely 
thought to be tidal ( |Dressler|1980l|Bicker et al.|2002|[Wilman et al.| 



In this paper we consider A 1367 which together with A 1656 
(the Coma cluster) comprises the Coma supercluster. A 1367 
(,z=0.022) lies at the intersection of two filaments; the first ex- 
tending 100 Mpc eastwards in the direction of the Virgo clus- 
ter an d the second running NE toward Coma (West & Blakeslee 
|2000| . A 1367, with a mean velocity « 6,240 km s~ and a ~ 
822 km s~ , is an unrelaxed Bautz-Morgan type II— III, spiral rich, 
6.9 x 10 14 M Q cluster (Boselli & Gavazzi 2006a) . Based on a red- 
shift to the cluster of 0.022 and Qm == 0.3 , Q,a = 0.7 and Ho = 
72 km s _1 Mpc -1 the angular scale is 1 arcmin = 24.8 kpc. 

A 1367 is, for several reasons, the perfect target to study en- 
vironmental effects on spirals. It has a much higher fraction of 
spirals than its more massive neighbour Coma. X-ray observa- 
tions of A 1367 carried out with ASCA, XMM-Newton (XMM) and 
Chandra indicate a dynamically young system with two principal 
sub-clusters (e.g., jDonnelly et al.|1998| . These two sub-clusters 
are themselves in the process of assembling from several smaller 
groups jCortese et al. 2004). This suggests that A 1367 may be 
more like the unrelaxed clusters at earlier epochs than the more 
relaxed systems typical of low redshift ( Oemler et al. 1997 ). Being 
a lower mass cluster of the kind in which [Poggianti et al. (2009) 
report enhanced rates of morphological transformation from spi- 
rals to SOs it is potentially a location in which to test the mecha- 
nism leading to galaxy transformation suggested by Wilm an~et al.| 
(2009). Furthermore, extensive optical and radio continuum stud- 
ies have been carried out for A 1367. Gavazzi, in a series of papers 
(Gavazzi et al. 2001a; [Gavazzi & Boselli[T999l|Gavazzi et al.|1995l 
|Gavazzi &"J affe 1987 ) reported several galaxies with radio contin- 
uum and Ha tails and exceptionally high star formation rates NW 
of the cluster centre, interpreted as arising from the interaction of 
galaxies with the cluster ICM. Also the recently discovered Blue 
Infalling Group (BIG), projected near the centre of A 1367 with 
a high relative velocity (Sakai et al. 2002; Gavazzi et al. 2003a; 
Cortese et al. 2006) , may constitute a unique laboratory for testing 
preprocessing theories in cluster assembly. 

This paper is the first stage of an investigation into the mech- 
anisms transforming late-type galaxies in Abell 1367 which is fo- 
cused on their H I content and position within the cluster. It is based 
on both ArecibcQsingle dish and NRACj^J-VLA synthesis H I imag- 
ing. Despite Arecibo's modest spatial resolution of ~3.5 arcmin, its 
homogeneous (in H I sensitivity) sky coverage, good sensitivity and 
high velocity resolution are well suited to study the global H I dis- 
tribution throughout nearby clusters. But interferometers, like the 
VLA, provide an order of magnitude better spatial resolution, and 
for nearby clusters the VLA's spatial resolution is sufficient to al- 
low the effects of environment on the H I distribution in individual 
galaxies to be studied. 

In this paper we present the results of a VLA low/medium 
resolution H I spectral line study of two adjacent fields in the central 
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sphere Centre, which is operated by Cornell University under a cooperative 
agreement with the National Science Foundation 

2 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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FIELD B 




2.1 VLA Field A - NW sub-cluster 



Figure 1. The positions of the two fields observed with the VLA are in- 
dicated by large circles showing the FWHM (32 arcmin) primary beams. 
VLA H I detections are shown with small black circles with crosses. The 
grayscale is a ROSAT X-ray image with white contours overlaid. The po- 
sitions of the SE, NW and W sub-clusters and the BIG compact group are 
labelled 



and NW parts of A 1367. The Hi imaging from these two fields 
covers a small fraction of the virial volume of the cluster. 

In section 2 we give details of the VLA H I observations in- 
cluding data reduction. Special attention is given to the method we 
developed to optimise the continuum subtraction, leading to im- 
proved signal-to-noise for low brightness H I features. We report 
our main observational results from the VLA H I imaging of indi- 
vidual galaxies in section 3 and, in section 4 the distribution of H I 
throughout the cluster volume from the AGES data. In section 5 we 
discuss the implications of our data for evolutionary mechanisms 
in the cluster, based on the H I content and colour distributions at 
both large and small scales as well as on the peculiar H I morpholo- 
gies shown by several galaxies. Our concluding remarks are given 
in section|6] 



2 OBSERVATIONS 

We used two sets of VLA H I observations both of which are well 



within the 104 arcmin (2.58 Mpc) virial radius of the cluster l |Moss| 
2006) and completely within the A 1367 AGES surveyed volume 
l Cortese et al.|2008) . Figure[T|indicates the FWHM primary beams 
for both VLA pointings (large circles). The figure also shows the 
intensity of the X-ray emission (ROSAT) from the cluster's ICM 
(white contours) and the positions of the VLA H I detections (small 
black circles enclosing crosses). Observational parameters for both 
fields are listed in Table[T] including the velocity resolution, central 
velocities, and the rms noise per field. 

The sensitivity is non-uniform between the two VLA fields be- 
cause of their differing array configurations and integration times. 
Additionally, for both fields, there is a sharp drop in sensitivity be- 
yond the primary beam FWHM radius, significantly impacting the 
three detections in Field A furthest from the pointing centre. 



Field A is centred near the NW sub-cluster (Cortese et al. 2004), 
-~10 arcmin NW of the cluster centre (position from NED between 
the NE and SE sub-clusters). Field A's position was chosen based 
on two considerations. First, the peculiar Zwicky infalling star- 
burst galaxies CGCG 097-073, CGCG 097-079, and CGCG 097- 
087 were imaged within the 32 arcmin FWHM of Field A's primary 
beam. Second, the first null of the beam was placed near 3C264 (a 
strong continuum source, with a flux density of 5 Jy at 1 .4 GHz) to 
minimise the effects of its sidelobes. 

This field was observed using the VLA in D configuration on 
the 26th and 30th of March 1999. Total observing time for the field 
was 15 hours. We used correlator mode 4 which gives two inde- 
pendent Intermediate Frequencies (IFs), each with dual polariza- 
tion (right and left-hand circular). The IFs for the first day were 
set so that each IF generated a position-position-velocity subcube 
(a, 8, velocity) with a velocity range of ~500 km s~ , with a small 
overlap between the velocity ranges. Combining the subcubes from 
both IFs produced a cube for the field with a velocity range of 
~1000 km s" 1 . The field was re-observed on the second day us- 
ing the same procedure but for a range of velocities adjacent to that 
observed on the first day. When combined, the subcubes from both 
days produced a single cube of 96 channels with a continuous ve- 
locity range from 5762 to 7810 km s _1 and a velocity width, using 
on-line Hanning smoothing, of 21 km s _1 per channel. The final 
cube's velocity range is more than twice the velocity dispersion of 
A 1367 (~822 kms^ 1 ). 

The data were calibrated and imaged following standard pro- 
cedures using the AIPS software package. For Field A, self- 
calibration with a single iteration was carried out to improve upon 
the standard complex gain calibration. This was necessary to mit- 
igate the effects of sidelobes from 3C264. We applied different 
robust weighting functions in the Fourier transform (AIPS task 
IMAGR), looking for a compromise between sensitivity and res- 
olution. Our final data cubes were produced with robust weight- 
ing (using ROBUST = 0) and have a resolution of about 45 arc- 
sec. The ROBUST option corrects the weights of the visibilities in 
the Fourier transform for the fact that there is a much higher den- 
sity of measured visibilities in the inner part of the wv-plane com- 
pared to the outer regions (Briggs 1995). This comes at a cost of a 
slightly increased noise compared to using what is known as 'nat- 
ural weights'. For this data set the Hi mass detection threshold is 
~7x 10 7 M© (corresponding to 3a over 2 consecutive 21 km sec -1 
channels). The equivalent column density sensitivity for emission 
filling the beam is then 1.9x 10 19 cm -2 . 

It was found that continuum subtraction was a critical step 
in the data reduction. Ordinarily one would search for line-free 
channels in each subcube and then use, e.g. UVLIN to subtract the 
continuum. The problem with our observations was that there was 
some H I emission at some position in almost every channel. This 
required a more elaborate approach consisting of deriving a con- 
tinuum map made up of those areas on each channel map with no 
line emission and averaging them. This map was subsequently sub- 
tracted from the line + continuum data to produce cubes containing 
only line emission. In practice this required the following proce- 
dure. First, continuum sources with peak flux density 10 mjy 
were modelled and removed in the wv-plane. This resulted in sub- 
cubes containing line and residual continuum sources of ^ 10 mJy. 
In the second stage AIPS tasks SQASH and COMB were used to cre- 
ate an average of all channels from the line plus residual continuum 
cube in the image plane, which was subtracted from the line plus 
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Table 1. VLA Observational Parameters 



Field 



Q'20(K) 



Array 
config 



Integration 
time 
[hours] 



Beam 
size 
["] 



Velocity 
resolution 
[kms" 1 ] 



Central 
velocity 
[kms" 1 ] 



[mJy beam 



mjy 



conversion 



[K] 



A 
B 



11 43 45.6 19 50 21 D 
11 44 50.0 19 47 00 C 



15.0 
3.6 



44x45 
15 x 15 



21 
11 



6800 
8200 



0.27 
0.65 



0.3 
2.7 



a Field A = VLA Project ID: AB900; Field B = VLA project ID AH801 



residual continuum cube. This enabled us to find line-free channels 
in the subcubes, although in most cases there were very few (just 1 
or 2). Continuum subtraction was then improved by repeating the 
subtraction but only using the average of these few line free chan- 
nels to produce subcubes with H I line emission only. The rms noise 
in these cubes was quite high because the continuum subtracted was 
based on very few channels. 

In the final stage these H I cubes were smoothed and blanked 
using AIPS tasks CONVL and BLANK, blanking out all areas in each 
channel containing line emission. These masks were in turn applied 
to the line + residual continuum subcubes to create, after again ap- 
plying SQASH, residual continuum-only maps for each of the sub- 
cubes. These residual continuum maps were then subtracted from 
the original line + residual continuum subcubes, resulting in sub- 
cubes with only H I emission, but with a much reduced rms noise. 
The continuum subtracted subcubes were subsequently combined 
with MCUBE. The final cube produced in this way had a noise of 
~0.27 mjybeam -1 which increased the number of Hi detections 
by 50 percent compared to the equivalent cube produced using 
UVLIN. 

2.2 VLA Field B - BIG 

The second set of VLA H I observations was centred on the Blue 
Infalling Group (Figure [TJ, overlapping slightly and located to the 
SE of Field A, with a central velocity of 8200 km s , i.e. the mean 
velocity of the BIG galaxies. Archival C-array observations (cor- 
relator mode 2 with on-line Hanning smoothing) from 8 December 
2002 were used to produce an image cube with 48 channels (af- 
ter discarding the noisy edge channels) with a velocity width of 
11 km s~ per channel. The velocity coverage was 7913-8587 
km s . Our final data cube has a resolution of ~15 arcsec. In 
the spatial overlap region between the two observed VLA fields, 
the velocity range extends from 5762 to 8587 km s _1 (except for a 
narrow gap at ~7900 km s ). For Field B, the equivalent H I mass 
detection threshold is ~8 x 10 7 Mq (corresponding to 3a in 2 con- 
secutive 1 1 km s _1 channels). This is equivalent to column density 
sensitivity for emission filling the beam of 2. 1 x 10 20 cm -2 . 



3 OBSERVATIONAL RESULTS 

3.1 The H I in individual galaxies: VLA observations 

Nine objects were detected in Field A (NW and W sub-clusters, 
|Cortese et al. 12004} and two more within Field B (BIG). Channel 
maps for each detected object are available in Appendix C online 
material. Figure [2] shows the VLA H I detections in the NW area 
of the cluster, and their locations relative to X-ray emission (from 
XMM-Newton). The XMM-Newton archive data were derived from 



an observation taken on 2001 Nov 22. After standard flare filter- 
ing, the exposure time was 20.6 ks for the pn and 29.9 ks for the 
two MOS cameras. Images in the 0.5-5.0 keV energy band were 
derived by combining the exposure-corrected images from all three 
cameras and convolving with a Gaussian. 

The morphological type, EW(Ho), VLA H I velocity includ- 
ing uncertainty, AVhi including uncertainty, H I mass, and offset 
(of the H I compared to the optical position) for the VLA detections 
are given in Table|2] 

The uncertainty in position for sources in the VLA fields is 
~ 4 arcsec (1.6 kpc) and 2 arcsec (0.8kpc) for Field A and B re- 
spectively, i.e, uncertainties are ~ l/10th of the synthesized beam. 
Within a velocity integrated intensity map of an individual spiral, 
positional uncertainly varies depending on S/N. This means posi- 
tional uncertainty within a map of a spiral is lowest in the region 
of velocity integrated intensity maximum. But we caution that the 
distribution beyond the H I intensity maximum region in spirals be- 
yond the FWHM of the primary beam is more uncertain. 

Four galaxies, CGCG 097-062, CGCG 097-068, CGCG 097- 
072 and CGCG 097-125, imaged with the VLA have their Hi in- 
tensity maximum offset relative to their optical counterpart (final 
column of Table [2} in directions which are not radial with respect 
to the cluster centre, i.e. these cases are inconsistent with a simple 
ram pressure scenario. The direction of these offsets may arise from 
a tidal interaction or a more complex ram pressure scenario such as 
gas fall-back (see subsection 3.3 for CGCG 097-072 below). 

Most of the galaxies have previously been detected in Hi 
with Arecibo jChincarini et al.|1983[ |Giovanelli & Haynes|1985[ 
|Gavazzi|1989||Springob et al.|2005|> a nd/or with the VLA ^Dickey] 
& Gavazzi 1991 ; Hota & Saikia 2007 ). None of the authors showed 
H I maps of their detected galaxies because of insufficient resolu- 
tion or sensitivity, except Hota & Saikia (20071 who presented Hi 
maps of the three most active starburst galaxies. 

All objects in Table [2] were also detected in AGES (Table [3] 
and |Cortese et al.|2008| l, with total fluxes and velocities generally 
in good agreement with ours. As an example we show, in Figure|3] a 
comparison of the AGES and VLA Field A spectra for a detection 
well within the VLA FWHM beam. Important exceptions to this 
agreement include CGCG 097-073 and A 1367 [GP82] 1292, which 
are confused in the AGES beam but are resolved by the VLA. 

While the determination of velocity for each channels is ac- 
curate to within a few km s _1 the range of velocities detected 
for each galaxy depends on sensitivity. The VLA detections of 
CGCG 097-068, CGCG 097-072 and CGCG 097-062 are beyond 
the Field A primary beam, where not only is the sensitivity lower 
than in AGES, but the fluxes are more uncertain due to the limited 
knowledge of the VLA primary beam beyond the half-power ra- 
dius. For example the noise at the position of CGCG 097-068 is ~5 
times greater than at the centre of the primary beam. 

In the following subsections we describe the most notable fea- 
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Table 2. Parameters for the VLA H I detected galaxies 



_X7I 


[7TI 


52000 




1771 

EW(Ha|^ 


WLK 


171 

<kvhi)|£| 


W 20 


av 2o M 
vix 


im 


171 

VLA 


m 

Offset[7j 




r h m s -| 
L J 


rol fl-, 
L J 




rAi 


[ km s _1 ] 


[ km s _1 ] 


[ km s -1 ] 


[ km s _1 ] 


[ km s _1 ] 


HO 9 Mr-il 




97-068 


11 42 24.5 


20 07 10 


Sbc 


41 


5958 


± 5 


346 


346 


±10 


6.9 


15" N 


GP1227 


11 43 13.0 


19 36 47 


Dw* 


30 


6239 


±11 


86 


86 


±22 


0.9 


none 


97-072 


11 42 45.2 


20 01 57 


Sa 


9 


6207 


± 4 


64 


64 


± 9 


0.4 


12" SE 


97-087 


11 43 49.1 


19 58 06 


Im 


77 


6738 


±27 


563 


649 


±54 


8.0 


30" NW 


97-079 


11 43 13.4 


20 00 17 


In- 


129 


7019 


±21 


216 


216 


±41 


1.3 


10" NW 


GP1292 


11 42 58.9 


19 56 12 


Dw* 




7203 


±23 


64 


64 


±46 


0.3 


none 


97-073 


11 42 56.4 


19 57 58 


SA:pec 


111 


7301 


±11 


173 


216 


±22 


2.0 


8" N 


97-091 


11 43 59.0 


20 04 37 


Sa 


23 


7377 


± 2 


259 


281 


± 4 


5.3 


none 


97-062 


11 42 14.8 


19 58 35 


Sa:pec 


37 


7723 


±13 


64 


64 


±26 


0.7 


15" SW 


97-125 


11 44 54.8 


19 46 35 


SOa 


23 


8158 


±18 


206 


206 


±35 


1.3 


12" SW 


K2 


1 1 44 50.6 


19 46 02 


Hll 




8158 


±12 


162 


162 


±24 







a Zwicky catalogue, except GP1277 = Abell 1367[GP82] 1227 and GP1292 = Abell 1367[GP82]1292 

^ Optical galaxy position from NED 

c Hubble classification from NED, * = our classification 

^ EW(Ha) is from GOLDMine, except GP1227 which is our measurement 

Velocity = the mean of the upper and lower velocit y defined by the line w idth (see note g) 



f Uncertainly for V H i = 1.5 (W2o-W 5 o)(S/N) -1 Schneider et al. 1990) , except where W20 = W50 in which case <x(V H i) 
= 1.5 (2 X channel width )(S/N) _1 . 

9 The velocity width based on the number of contiguous channels with ^2 a detections. 
^ Uncertainty in the line width = 2<t(Vhi) 

1 M H i = 2.36 x 10 2 D 2 S H i where M H i is in M0, D in Mpc and Shi is the VLA flux in mjy km s 1 . Mass based on a 
distance of 91.6 Mpc (V A 1367 = 6595 km s" 1 from NED and H Q = 72 km s" 1 Mpc -1 .) 

3 The offset direction is the projected position of the Hi intensity maximum relative to its optical counterpart, except for 
CGCG 097-087 where the direction and distance are taken from the NIR counterpart. The projected distance between the two 
intensity maxima is given in arcsec with a positional uncertainty of ru 4 arcsec. 
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Figure 2. NW region of A 1367: X-ray emission (XMM) image with white 
contours. Black contours indicate H I emission from the VLA H I detections 
in field A, with the outer contour indicating a column density of N hi = 3 x 
10 19 cm -2 . Values in parentheses below the Zwicky identifier are optical 
velocities in kms -1 . 




Figure 3. CGCG 097-087 VLA (solid line) and AGES (dashed line), inte- 
grated H I spectra. The feature at 7000 km s~ 1 in the AGES spectrum is an 
RFI artifact. 



tures of the VLA detections. The images displayed here were pro- 
duced after applying a primary beam correction. The a and S axes 
of all maps are in J2000 coordinates. 



3.2 The star-forming galaxies (CGCG 097-087, 
CGCG 097-073 and CGCG 097-079) 



These three starburst galaxies have been extensively studied in the 



past (Gavazzi et al. 2001a,b; jNulsen 


1982; Gavazzi et al.IT984; 


Gavazzi & Jaffe|1987| Gavazzi 


1989 


Gavazzi et al.|1995||Dickey| 


& Gavazzi|1991||Boselli et al. 


1994 


Hota & Saikia|2007|(. They 



all have exceptionally high star formation rates, EW(Hq + [NII])> 
75A (Table bl and display spectacular radio continuum and Ha 
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tails pointing away from the cluster centre (Gavazzi et al. 2001b; 
|Gavazzi &"j affe 1987) . Some observational features of these three 
objects have been interpreted as effects of ram pressure stripping 
produced by the ICM, in particular the enhanced star formation 
activity and the displacement of atomic hydrogen relative to their 
stellar disks, although in the case of CGCG 097-087 much of the 
enhanced star formation rate is probably attributable to a merger 
(Gavazzi et al. 2001a; |Martig & Bournaud||2008| > The Hi in the 
central parts of these starbursts has been recently mapped by |Hota| 
|& Saikia] < |2007) ). The maps presented here, though, show additional 
large scale H I features as a result of our improved continuum sub- 
traction (see Sect. |2. 1 [ l. 

• CGCG 097-087: The integrated Hi map (Figure [4} shows 
the Hi intensity maximum is offset by ~30 arcsec (12 kpc) to 
the NW of the nucleus, confirming the H I is asymmetrically dis- 
tributed, as earlier reported by |Gavazzi] (19891 and |Dickey &| 
|Gavazzi| ( |1991} . Higher spatial resultion and lower sensitivity Hi 
maps have previously been presented for the inner parts of this 
galaxy by |Hota & SaiMa| ( (2007) l. The main body of the galaxy dis- 
plays rotation in the velocity range ~6500-6900 km s _1 (Fig.[5]l. 
We do not see the jump in velocities observed in Ha and CO 
(Gavazzi et al. 2001a; |Boselli et aT]|1994) , but this may be due 
to our low spatial resolution. We are not able to confirm the merger 
hypothesis suggested by the Ha kinematics (Gavazzi et al. 2001a 
et al. ; |Amram et al.|20 02), but the H I velocity field has a much 
steeper velocity gradient SE of the nucleus, coinciding with the 
disturbed optical disk. We also see that the H I disk is sharply trun- 
cated, with its SE edge approximately coinciding with the edge of 
the optical disk. An impressive H I tail, coinciding with the radio 
continuum tail ( Gavaz zTet al.|1995) , extends at least 70 kpc to the 
NW of the nucleus, with an increasing offset from the plane of the 
optical galaxy. At the end of this tail the H I map shows a clump 
without an optical counterpart. The clump has a velocity width of 
~100kms _1 . 

• CGCG 097-073 and A 1367 [GP82]1292: Figure § shows 
the Hi in CGCG 097-073 having a slightly asymmetric distribu- 
tion, the intensity maximum being displaced ~8 arcsec (3 kpc) to 
the north of the optical disk, confirming the earlier report of H I 
asymmetry (Dickey & Gavazzi 1991 1. The Hi displacement is in 
the same direction as the Ha and radio continuum tails (Gavazzi et 
al. 1995, 2001b). H I does not appear truncated, which is consistent 
with the normal Hi content (Hi deficiency = 0.02). Also shown 
in the same figure is an Hi detection at the position of the blue 
dwarf galaxy Abell [GP82]1292 which appears undisturbed. Given 
the projected distance of ~ 1.9 arcmin (47 kpc) and the 1 1 km s -1 
velocity separation between the two (perturbatioij^] parameter p gg 
~ 0.004), a gravitational interaction cannot be ruled out. 

• CGCG 097-079: H I has been previously detected in this 
galaxy (Gavazzi 1989, Dickey & Gavazzi 1991). Figure [7] shows 
the H I to be symmetrically distributed but with the intensity max- 
imum displaced by ~10 arcsec (4 kpc) to the NW of the optical 
disk, as previously reported by Gavazzi ( 1989) and H ota & Saikia| 
(2007). The H I offset is in the same direction as the Ha and radio 
continuum tails (Gavazzi et al. 1995, 2001b) . 



3 Pgg = ^ M """ p ^ A ^ a '^ where M„ a ; and M comp are the masses of the 
galaxy and companion respectively, d is the separation and r is the galaxy 
disk radius jByrd & Valtonen|1990} . 
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Figure 4. CGCG 097-087: White contours are from a robust H I surface 
density map, with the outer contour indicating a column density of Njji 
= 3 x 10 19 cm -2 , with higher levels at 8, 20, 40, and 80 x 10 19 cm -2 , 
overlaid on an SDSS ('-band image. The first contour also corresponds to a 
8 a detection in three channels. The size of the D-array beam is indicated 
with the white circle. 



+ 19°59 



+ 19°57' 




ll n 43 m 56 : 



Figure 5. CGCG 097-087 velocity field: White contours show velocities 
from 6500 to 7000 km s , in 50 km s _1 steps based on the robust 
cube overlaid on an SDSS r'-band image. The size of the D-array beam is 
indicated with the white circle. 



3.3 Giant late-type galaxies (CGCG 097-062, 

CGCG 097-068, CGCG 097-072 and CGCG 097-091) 

In addition to the three starbursts we were able to map four more 
spirals in the NW and W sub-clusters (Figure[2](. 

• CGCG 097-062: This optically asymmetric Sbc (NED) or 
Pec Galaxy OnLine Database Milano Network (GOLDMin^J, 



http://goldmine.mib. infn . it/ 
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Figure 6. CGCG 097-073 (upper contour) and A 1367 [GP82]1292 
(lower contour): Contours trace the robust H I surface density, with the 
outer contour indicating a column density of Nhi = 3 x 10 19 cm -2 , with 
higher levels at 10, 17 and 24 x 10 19 cm" 2 overlaid on an SDSS /-band im- 
age. The first contour also corresponds to a 7 cr detection in three channels. 
The size of the D-array beam is indicated with the white circle. 




Figure 7. CGCG 097-079: White contours are from a robust Hi surface 
density map, with the outer contour indicating a column density of Nhi 
= 3 x 10 19 atoms cm -2 , followed by contours at 10 and 17 x 10 19 atoms 
cm -2 on an SDSS ('-band image. The first contour also corresponds to a 
7 cr detection in three channels. The size of the D-array beam is indicated 
with the white circle. 
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Figure 8. CGCG 097-062: White contours are robust H I surface density 
levels with the outer contour corresponding to a column density of Nhi = 3 
x 10 19 atoms cm -2 , with higher levels at 8 and 13 x 10 19 atoms cm~ 2 , on 
an SDSS !-band image. The first contour also corresponds to a 2 cr detection 
in three channels. The size of the D-array beam is indicated with the white 
circle. 



galaxy is well beyond Field A's FWHP beam, but our weak Hi 
detection is consistent with the Hi deficiency (AGES) of 0.35. Fig- 
ure [8] shows the Hi intensity maximum displaced ~15 arcsec (6 
kpc) to the SW of its optical counterpart, in the direction of the 
optical tail. CGCG 097-062, like CGCG 097-087 has a pronounced 
asymmetric H I spectrum (Gavaz zi|l989| >. 

• CGCG 097-068: Despite its position well beyond the 
FWHP region of the VLA's primary beam, CGCG097-068's Hi 
signal was sufficient to show a double-peaked VLA spectrum 
which approximately matches its AGES counterpart, strongly sug- 
gesting that this massive Sbc galaxy is inclined and Hi rich (Hi 
deficiency of -0.28, Table BJ. This galaxy was tentatively reported 
as having its H I intensity maximum displaced north of the optical 
disk {Dickey & G avazzi 1991 1. The position of the VLA H I inten- 
sity maximum ~15 arcsec (6 kpc) in projection north of the optical 
nucleus confirms this (Figure[9j. 

• CGCG 097-072: This Sa galaxy lies, in projection, close to 
the W sub-cluster galaxy density maximum ( [Cortese et al.| 2004l. 
Our Hi detection is weak, consistent with its high Hi deficiency 
of 0.55 (AGES). Figure[T0| shows an Hi intensity maximum offset 
~12 arcsec (5 kpc) to the SE of the optical nucleus, and suggests 
a truncated Hi disk. However, this requires further confirmation 
as the galaxy is beyond the FWHP beam of Field A. The high Hi 
deficiency and anomalous direction of its H I displacement relative 
to that seen in nearby spirals is similar to NGC 4848 in Coma where 
a truncated H I disk was reported by Bravo- Alfaro et al. ( 2001} and 
further CO-imaging and numerical simulations by Vollmer et al. 
(2001b) revealed a case of gas fallback. 

• CGCG 097-091: Hi has been previously detected in this 
galaxy with Arecibo and the VLA ( [Dickey &~G avazzi 199T}, but 
was not spatially resolved. The AGES Hi deficiency of -0.23 shows 
the galaxy to be Hi rich. CGCG 097-091 has the most symmetric 
H I morphology of our VLA detections, with a near perfect coin- 
cidence of position and velocity between the optical nucleus and 
the Hi intensity maximum, (Tables [3] and [2} . These characteristics 
together with the velocity field (FigurefTT), which shows normal ro- 
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Figure 9. CGCG 097-068: White contours are from a robust H I surface 
density map, with the outer contour indicating a column density of Nhi 
= 10 x 10 19 atoms cm -2 , higher levels at 30, 50 and 70 x 10 19 atoms 
cm -2 on an SDSS f'-band image. The first contour also corresponds to a 5 
a detection in three channels. The size of the D-array beam is indicated 
with the white circle. 
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Figure 11. CGCG 097-091 H I velocity field: based on a robust cube, 
Contours are 7275, 7300, 7325, 7350, 7375, 7400, 7425 and 7450 km s" 1 
on an SDSS r— band image. The size of the D-array beam is indicated with 
the white circle. 
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Figure 10. CGCG 097-072: White contours are from a robust H I surface 
density map, with the outer contour indicating a column density of Nhi = 

3 x 10 19 atoms cm -2 and the higher contours are 5 and 7 x 10 19 atoms 
cm -2 , on an SDSS /-band image. The first contour also corresponds to a 

4 a detection in three channels. The size of the D-array beam is indicated 
with the white circle. 



maximum is offset ~12 arcsec (5 kpc) westward from the optical 
centre (Figure[T2j. The high column density H I to which Field B is 
sensitive is asymmetrically distributed with its major axis running 
approximately SE - NW. Comparing our H I images to the earlier 
more sensitivity observations by Sakai et al. (20021 with the West- 
erbork Synthesis Radio Telescope (WSRT) indicates that BIG con- 
tains extensive diffuse H I, which our observation was not sensitive 
enough to detect. Comparison between our H I spectrum and AGES 
confirms the presence of this diffuse H I in the velocity range from 
8000 to 8500 km s _1 . Further details of our observation of BIG 
can be found in Appendix A. 

3.5 Dwarf galaxy A 1367 [GP82] 1227 

This isolated irregular dwarf galaxy is detected in H I, both by the 
VLA and in AGES, to the SW of the cluster centre in a region oth- 
erwise devoid of Hi detections (Figure [1"3|>. The VLA and AGES 
H I spectra are in good agreement and display a single narrow peak 
(W 50 =56 ±8 km s~ x ; Cortese et al. 2008), typical of a dwarf ir- 
regular galaxy. The VLA observations show the H I intensity max- 
imum is offset 3 arcsec (1 kpc) to the south of its optical counter- 
part. Details of a low-resolution optical spectrum obtained for this 
galaxy can be found in Appendix B. 



tation, combine to suggest CGCG 097-9 1 is a normal spiral which 
is not interacting with its environment in any significant way. 



3.4 The H I around the Blue In-falling Group (BIG) 

Our VLA H I observation of the BIG (Field B) has a higher spa- 
tial resolution (~15 arcsec ), compared to Field A (~45 arcsec), 
but at significantly lower sensitivity because of the shorter integra- 
tion time. The strongest detection in Field B is near the disturbed 
SOa galaxy CGCG 097- 125 and we confirm that the Hi intensity 
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Figure 12. CGCG 097-125 and knot near K2: White contours are from a 
robust H I surface density map, with the outer contour indicating a column 
density of Nhi = 4 x 10 20 cm -2 , higher levels are at 8 and 12 x 10 20 
cm -2 on an SDSS g-band image. The first contour also corresponds to a 
4 a detection in three channels. The size of the C-array beam is indicated 
with the white circle. 
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at 10 and 17 x 10 cm" , overlaid on an Ha image. The first contour also 
corresponds to a 4 a detection in three channels. The size of the D-array 
beam is indicated with the white circle. 



4 H I (AGES) AND COLOUR DISTRIBUTION IN A 1367 

To place our VLA observations in context we used the AGES H I 
survey together with SDSS data to analyse the Hi content and 
colour of bright late-type galaxies throughout the central volume 
of the cluster (including the VLA fields). 

As part of the AGES project the 305-m Arecibo Telescope was 
used to survey Hi in a volume centred on A 1367 covering about 

5 square degrees, between ll' l 34 m 00 s < a 2 ooo < ll fe 54 m 15 s , 
19 d 15 m < 5 2 ooo <20 d 20 m , in a velocity range ~1 100-19000 



method from 



Figure 13. [GP82]1227:White contours are from a robust Hi surface 

density map. The outer contour is N jj j = 3 x 10 19 cm -2 , higher levels are ISolanes et al I 

.» in -_J n_ i/»19 -2 i-j it. : tl. . I — 



km s _1 i Cortese et al.pOOl i. The survey has an angular resolu- 
tion of 3.3 x 3.8 arcmin and a velocity resolution of ~10 km s -1 
with a mass sensitivity limit of 6 x 10 8 M© of H I over a velocity 
W50 = 200 km s _1 . Positional uncertainties for the detections are 
estimated to be ~ 18 arcsec (~7 kpc). The observations were made 
in drift scan mode providing a consistent flux limited H I detection 
threshold throughout the AGES volume. The observations and data 
reduction process is descr ibed in|Auld et al.H2006) . More details 
on this survey are given in Cortese et al. | 2 008| >. 

In the field, significant amounts of Hi are expected to be 
present in the disks of both spiral and dwarf galaxies. In the case of 
spirals the amount of H I can be approximately related to the opti- 
cal disk diameter and the galaxy's morphological type ( Ha ynes &] 
|Giovanelli| 1984) |Solanes et al.|1 996), but for dwarfs the H I mass 
fraction is both larger and more uncertain. The AGES detection 
limit implies that spiral galaxies in A 1367 with an optical disk di- 
ameter < 0.37 arcmin (~9 kpc) would not be detected in AGES 
( |Solanes et al.|1996| |Cortese et al.|2008) . Galaxies with disk di- 
ameters below this limit are almost certainly dwarfs, so that AGES 
is essentially only sensitive to Hi in late-type spiral galaxies of 
A 1367. However two dwarf galaxies were detected within our sam- 
pled volume, showing that AGES includes galaxies from the high 
end of the dwarf H I mass function. 

AGES allows us to investigate the Hi deficiency of A 1367's 
spirals, down to its sensitivity limit, by comparing their expected 
H I content with actual AGES determined H I content. To do this a 
catalogue of A 1367's brightest spiral members was compiled, with 
the following criteria. We intially selected late-type galaxies with 
an SDSS <?-band magnitude < 15.5 within the sampled volume, 
i.e., 1 square degree surrounding the NW sub-cluster core in the 
velocity range 4000-9000 km s , (approximately six times the 
velocity dispersion of the cluster). Only those objects with an op- 
tical disk diameter > 0.37 arcsec from GOLDMine (Gavazzi et al. 
2003b) or, if unavailable, from the NASA Extragalactic Database 
(NED), were selected. For the selected late-type galaxies we de- 
termined their Hi deficiency by comparing expected and observed 
H I content using the method of Haynes & Giovanelli ( 1 984 ). Hub- 
ble type was taken from GOLDMine or, if unavailable, from NED. 
We selected all galaxies with Hubble type later than SO. The only 
exception, CGCG 097- 125, an SOa, was included because of its un- 
usually large H I mass. 

Applying these criteria we found the 26 bright late-type galax- 
ies listed in Table[3]which gives details of their SDSS g band mag- 
nitude and g-i colour, optical diameter from GOLDMine/NED, 
Hi-mass based on AGES flux, Hi deficiency (or lower limit for 
AGES non-detections) based on AGES flux calculated using the 
Haynes & Giovanelli (1984) and parameters from 



1996} , and separation parameters for galaxies with a 
close companion. A g-i colour of 1 . 1 approximately coincides with 
the red sequence threshold displayed in the colour-magnitude plot 
in Figure 12 of Cortese et al. (2008 ). For our sample, those galaxies 
with g-i colour < 1.1 are referred to throughout the paper as blue 
galaxies and the rest referred to as red galaxies. 

Figure [14] shows the optical positions of the 26 spirals (red 
and blue symbols) and the positions of the 18 AGES H I detections 
(white circles). In the figure the positions of the 16 selected galaxies 
with blue SDSS g-i colours are indicated with blue symbols and 
the 10 with red SDSS g-i colours with red symbols. Figure [T4] also 
indicates the distribution of ICM gas in the cluster based on ROSAT 
X-ray observations (black contours). 

A shift in position in Fig. [14] between an AGES H I detec- 
tion (white circle) relative to its respective optical position (red and 
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Figure 14. Positions of selected late-type galaxies in the volume sampled. The colour of the symbol (red or blue) used for each spiral indicates its SDSS g—i 
colour: blue (Sjl.l), red (>1.1). AGES H I detections are marked (white circles). The symbol used for each spiral indicates its evolutionary state (discussed in 
section [53) as follows; A-(asterisk), B-(open circle), C-(filled circle) and D-(open square) BIG spirals-(cross) and unclassified-(diamond). 

ROSAT X-ray intensities (black contours) are overlaid on an SDSS ('-band image. Details for the galaxy identifiers are given in note (a) of Table[3] 



blue symbols) is an indication of a real displacement when the off- 
set exceeds 18 arcsec (7 kpc), the uncertainty in the AGES posi- 
tions. This is only the case for CGCG 097-087 and CGCG 127- 
032. However the much higher pointing accuracy of the VLA, 
~4 arcsec (1.6 kpc), reveals that CGCG 097-062, CGCG 097-068, 
CGCG 097-072, CGCG 097-073, CGCG 097-079, CGCG 097-087 
and CGCG 097- 125 all have Hi intensity maxima offsets relative 
to their optical counterparts larger than 3 a. In all cases the offset 
direction is consistent with that from the higher uncertainty AGES 
data. The projected magnitude and direction of the VLA offsets 



given in Table|2]is consistent with that expected for spirals on radial 
orbits experiencing ram pressure stripping, except CGCG 097-062, 
CGCG 097-068, CGCG 097-072 where the offset direction is not 
radially with respect to the cluster centre. 



The four AGES detections in Figure [14] without a selected 
late-type galaxy counterpart are the SO galaxy CGCG 127-032, 
the two dwarfs, [GP82] 1227 and [B085] 079, and an AGES 
detection to the NE of CGCG 097-068. CGCG 097-1 19 was not 
detected in the AGES survey because of its proximity to 3C264 
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Table 3. Parameters for the brightest late-type galaxies in A 1367 
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14.92 


1.14 


0.45 


Irr 


^0.21 


194 24 


C 


97-121 


11 44 47.04 


20 07 30.32 


6571 


12.86 


1.18 


1.20 


Sab 


0.38 




C 


97-114 


11 44 47.8 


19 46 24.31 


8293 


15.12 


0.57 


0.54 


Pec 


^0.34 


22 102 




97-119 


1 1 44 47.95 


19 41 18.56 


4895 


14.18 


0.98 


0.60 


Sa 


^0.37 






97-120 


11 44 49.16 


19 47 42.14 


5595 


12.53 


1.12 


1.32 


Sa 


^0.78 




D 


97-122 


11 44 52.23 


19 27 15.12 


5468 


13.65 


0.93 


1.45 


Pec 


0.47 




B 


97-125 


1144 54.85 


19 46 35.18 


8271 


13.81 


1.20 


0.84 


SOa 


>0.14 


22 102 




97-129-NED1 


11 45 03.88 


19 58 25.24 


5085 


12.29 


1.21 


2.36 


Sb 


-0.07 




C 


97-138 


1145 44.6 


20 01 50.8 


5313 


13.99 


0.42 


0.75 


Pec 


-0.12 




A 


B085 074 


11 45 06.57 


20 09 23.21 


6151 


14.69 


1.03 


0.38 


S 


0.08 




B 


97-129 NED2 


11 45 06.95 


19 58 0.66 


7546 


14.66 


1.12 


1.07 


Sbc 


^0.8 




D 



Zwicky catalogue, except [BO85]074 from Butcher & Oemler ( 1985 , with the suffixes indicating the NED pair identifier. 
Optical galaxy positions ex SDSS DR6 
from NED 
ex SDSS DR6 
ex SDSS DR6 

ex GOLDMine except CGCG 097- 1 1 1 and [BO85]074 which are taken from NED 
ex GOLDMine except [BO85]074 which is ta ken from NED 

and the parameters from 



Hi deficiency is calculated using the method of Haynes & Giovanelli 



1984 



Solanes et al. 



1996 



using total H I flux measurements from AGES and galaxy diameters equal to the major axis diameter from GOLDMine 
(except [BO85]074 which is taken from NED). For 97-079 the mass was calculated from our VLA total flux because the 
AGES flux was affected by RFI. For AGES non-detections, the Hi deficiency lower limit was calculated using the AGES 
detection limit of 6 x 10 s Mq. No Hi deficiency could be calculated for CGCG 097- 119 because flux from this galaxy was 
confused in the AGES beam with 3C264. The calculated H I deficiencies are highly uncertain because of the intrinsic spread 
in H I properties of individual spirals and we adopt an uncertainty of ± 0.3 dex throughout this paper. 

* A galaxy is considered to be a likely companion if its separation is 2 arcmin, the velocity difference is ^ 500 km s — 1 
and the companion has a diameter ^0.5 a rcmi n. 
3 Evolutionary state described in section .' 



5.3 



and is therefore not included in our analysis. 



ICM and the evolutionary state of the spirals. 



5 DISCUSSION 

The projected distribution, colour, H I content and velocities of the 
selected spirals throughout the sampled volume on large and small 
scales provides information on the recent history of the cluster as 
well as clues to the mechanisms impacting the galaxies' ISM. In 
the following subsections we consider the global Hi properties, 
the role of ram pressure stripping attributable to the NW subcluster 



5.1 Global H I properties 

Considered on a scale of several megaparsecs the pattern of H I 
deficiencies projected in RA for Coma and A 1367 in Figure 2 of 
|Gavazzi| (JT989 1 show striking differences. From their figure it is 
clear that with few exceptions the Coma spirals across the central 
~5° of cluster are significantly H I deficient (> 0.3). In contrast for 
A 1367 the significantly Hi deficient spirals are found in a more 
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Figure 15. The Hi deficiency of all the selected late-type galaxies in A 
1367, (excluding CGCG 097-119) as a function of projected distance from 
the cluster centre. The following symbols are used to indicate the state (dis- 
cussed in section |5~3) of each spiral: A-(asterisk), B-(open circle), C-(filled 
circle) and D-(open square). For AGES non-detections the Hi deficiency is 
a lower limit and is indicated with arrowheads. The two BIG spirals are con- 
fused in the AGES beam (triangles with upper error bar only) and in both 
cases the value is the lower limit for Hi deficiency. The region between the 
horizontal lines approximates the natural range of variation expected in the 
field population. 



narrow range of projected RA (~1°) and even at its centre a major- 
ity of the spirals are not deficient (H I deficiencies < 0.3) 

Within our sampled volume for A 1367 (~ central 1.5 
Mpc) the bright spirals, irrespective of Hi content, are found 
preferentially projected onto the northern half (Figure [14} , with the 
AGES detections (white circles in Figure [14) highly concentrated 
in a band running roughly E-W between declinations of 19°55' to 
20° 15'. 

As an initial step to understand the global H I properties in the 
sampled volume we plotted, in Figure^] the Hi deficiency of all 
the selected late-type galaxies as a function of projected distance 
from the cluster centre, with lower limits for H I non-detections, 
and BIG members confused in the AGES beam. The plot shows 
that the mean Hi deficiency of the selected galaxies is appreciably 
greater than in the field indicating that the cluster environment has 
impacted the ISM of a significant fraction of them. However, un- 
like Coma where a clear increase in Hi deficiency is seen within a 
radius of ~1 Mpc of the core ( |Bravo-Alfaro et aLp OOO), there is no 
obvious correlation between Hi deficiency and projected proximity 
to the cluster core. It should be noted, though, that the two inter- 
acting major sub-clusters revealed by the X-ray observations (see 
Fig.[T]( complicate the investigation of any correlation between the 
cluster's ICM density and Hi deficiency and in section [5T2| we try 
to address the effects which are specific to the NW sub-cluster. 



5.2 Ram pressure stripping by the NW subcluster ICM 

We are interested to know if the observed H I deficiencies can be 
explained in terms of ram pressure stripping. In order to quantify 
the role of this mechanism around the NW sub-cluster, we car- 




dm, (Mpc) 

Figure 16. Ram pressure as a function of the distance from the X-ray emis- 
sion maximum of the NW sub-cluster for v re i of 500, 1000, 1500 and 
2000 km s _1 . The dashed line indicates the threshold for ram pressure 
stripping based on Cayat teet al.H1994) . 



ried out a simple estimate of the strength of ram pressure as a 
function of radius from the NW sub-cluster ICM core. The ICM 
density (picm) distribution was estimated using the hydrostatic- 
isothermal /3-model of Cavaliere & Fusco-Femiano (1976| with 



the gas profile parameters from Donnelly et al. ( 19981. Ram pres 



sure values were derived from Gunn & Gott's (1972) equation: 



= PICM X V r 



where picM is the density of the ICM 



at the galaxy position and v re i is the galaxy's velocity relative to 
the ICM. We have no information about Vrei for the selected galax- 
ies, although we consider a galaxy's radial velocity, relative to the 
cluster systemic velocity, as a lower limit for v re i. 

The results of this modelling of P ram as a function of prox- 
imity to the NW core for four different values of v re i, up to 
2000 km s _1 (approximately the radial velocity difference of the 
BIG), are shown in Figure [T6] A spiral near the NW sub-cluster 
is expected to suffer ram pressure stripping if P Tam values exceed 
the restoring force of the galaxy's stellar and gas disks, at the gas 



disk's periphery, defined as F r = 27r G a stellar a g 



The thresh- 



old ram pressure required to remove Hi from the outer edge of 
the disk can be estimated from the typical value for the restoring 
force per unit area of F r — 2 x 10~ 12 dyncm -2 , obtained by 



|Cayatte et al.| {1994) for spiral galaxies in Virgo (dashed line in 
Figure 1 16). The crossing time for the cluster is 1.7 Gyr (Boselli & 
Gavazzii 2006a) making it likely that spirals within the central Mpc 
have been subject to ram pressure over periods of several times 10 8 
yrs. Modelling by Roediger & Briiggen (2007) indicates that for 
a mid-size spiral (V ro t = 200 km s~ ) sustained exposure to ram 
pressures over periods of the order of 10 s yrs of 10~ 12 dyn cm -2 
will result in only minor stripping of Hi, pressures of 10 _1 dyn 
cm -2 will cause significant stripping, while 10" 10 dyn cm -2 will 
almost completely remove the Hi. Massive galaxies are expected 
to suffer proportionately less H I removal, and less massive galaxies 
proportionately more. 

The implication of Figure [T6] and the Roediger & Briiggen 

(2007 i models is that any galaxy with v re i ^ 1000 km s~ any- 



where within a radius of 1 Mpc of the NW sub-cluster will undergo 
weak to moderate ram pressure stripping. But removal of large frac- 
tions ofHl, i.e. Pram 1? 10 _ 1 1 dyn cm~ 2 , would require a transit 
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Figure 17. The Hi deficiency of selected state A, B and C galaxies as a func- 
tion of distance from the NW-ICM core. The following symbols are used 
to indicate the state (discussed in section [53) of each spiral: A-(asterisk), 
B-(open circle) and C-(filled circle). 

Galaxies with only upper error bars (arrowheads) are AGES non-detections 
so the value is the lower limit for Hi deficiency. The region between the 
horizontal lines approximates the range expected in the field population. 
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Figure 18. The Hi deficiency of all the selected late-type galaxies (exclud- 
ing CGCG 097-1 19 and BIG spirals) plotted against their SDSS g-i colour. 
The long dashed lines mark the boundaries of the parameters defining the 
four evolutionary states marked A to D. The following symbols are used 
to indicate the state of each spiral: A-(asterisk), B-(open circle), C-(filled 
circle) and D-(open square). For AGES non-detections the Hi deficiency 
is a lower limit and is indicated with arrowheads. The region between the 
dotted horizontal lines approximates the natural range of variation expected 
in the field population. 



of the high density ICM within ~ 0.75 Mpc of the cluster centre at 
a velocity in excess of 1500 km s _1 or transit within ~ 0.25 Mpc 
of the core at 1000 km s - 1 . 

To test the model predictions we selected the 14 spirals from 
our sample which were potentially suffering ram pressure stripping 
solely attributable to the ICM of the NW sub-cluster. Virtually all 
of these galaxies have radial velocities relative to the cluster of 
>500 km s _1 and it seems reasonable to assume that most, if 
not all, are suffering ram pressure stripping of varying strengths 
given they are at projected radii of ^ 0.8 Mpc. Figure [IT] shows a 
plot of Hi deficiency against distance for these galaxies and, in- 
terestingly, only a modest increase in Hi deficiency is seen when 
compared with the field. Furthermore, no correlation is seen with 
proximity to the NW sub-cluster ICM core. This observational re- 
sult suggests that ram pressure stripping associated with the NW 
sub-cluster ICM core is relatively weak. The presence of undis- 
turbed H I in CGCG 097-073 's dwarf companion projected close to 
the NW-ICM core provides additional support to this statement. 

Although the limited evidence for ram pressure stripping is 
broadly consistent with the expected ram pressure at the projected 
distances to the the NW sub-cluster core (Figure [T6), the lack of 
correlation with radial distance from the core may imply that galax- 
ies seen in projection to be located close to the core are in reality 
spread over a considerable distance. In the case of CGCG 097-068, 
a tidal interaction cannot be ruled out as an explanation of the ob- 
served H I characteristics. 



5.3 Evolutionary states of the spirals 

The colour and Hi deficiency of a spiral are both properties in- 
dicative of its evolutionary state. In Figure [T8] we distinguish four 
broad evolutionary states (A-D). We classify the highly Hi defi- 
cient galaxies (Hi deficiency lower limit of JsO.7 ) as state D spi- 
rals. The spirals with moderate or no Hi deficiency (i.e., Hi defi- 



ciency ^0.7) are classified into three states (A, B and C) based on 
their SDSS g-i colour. C state spirals are red as defined in section 
[4] (i.e., g — i ^ 1.1). The remaining blue spirals with moderate or 
no Hi deficiency were classified as state A if their SDSS g-i colour 
was <0.76 (the colour of the bluest type of late-type field spirals 
from |Durbala et a l. (2008 ), their Table 3) and state B for those with 
0.76 ^g-i <1.1. Table [4] summarises the parameters used to deter- 
mine the evolutionary state of each spiral and the symbols used for 
them throughout the paper. The evolutionary state determined for 
each spiral is shown for convenience in the final column of Table|3] 
CGCG 097-125 and CGCG 097-1 14 have been excluded from the 
analysis because their evolution is likely to have been dominated 
by membership of the BIG compact group rather than the cluster. 

While the choice of boundaries of our states could be debated, 
it is clear that the cluster contains spirals at significantly different 
stages of evolution, with the spirals in each state appearing to have 
reached a broadly similar evolutionary stage, probably as a result 
of sharing or passing through a common environment. 

The discussion below is focused on the more evolved State C 
and D spirals where the impacts of the cluster environment are more 
pronounced. In contrast spirals with less evolved states A and B 
are probably at earlier stages of their interactions and consequently 
show milder impacts, although a key question remains to be ad- 
dressed, i.e. the mechanism(s) triggering the starbusts (state A). 

The H I deficiencies of AGES detected state C spirals are mod- 
est compared to lower limits for the state D spirals, suggesting that 
they have not passed through the cluster's highest density ICM. 
The state C spirals include 2 of the 3 potentially interacting pairs, 
where the companion is projected within 2 arcmin (~50 kpc), has 
a velocity within 500 km s _1 and both galaxies are ^ 0.5 ar- 
cmin (~ 12 kpc) in diameter (Table[3j. The perturbation parameters, 
p 99 , for the state C pairs are ~ 8.7 (CGCG 097-102), and ~ 1.8 
(CGCG 097- 111) significantly above the value where the ISM is 
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Table 4. Evolutionary states of spirals 



Plot H I g-i colour 

State Hi Colour symbol number def 

deficiency SDSS g-i galaxies median median 



A 


<0.7 


<0.76 


* 


6 


0.14 


0.59 


B 


<0.7 


0.76 <and <1.1 


o 


7 


0.47 


0.99 


C 


<0.7 




• 


4 


0.30 


1.20 


D 


>0.7 


any 


□ 


6 


0.79 


1.15 



expected to be impacted. The interaction time scales for the pairs is 
of the order of 10 8 years, in agreement with the value from Boselli 
& Gavazzi 2006a. It is important to note that Hi emission from 
close pairs could potentially be confused within Arecibo's ~3.5 
arcmin beam underestimating the Hi deficiency (AGES) for each 
member of the pair. But in both cases the spiral's companion is an 
early type making it unlikely the companion's H I makes a signifi- 
cant contribution. Moreover three state C spirals (CGCG097-102, 
CGCG 097-1 1 1 and CGCG 097-121) have optically disturbed mor- 
phologies in SDSS i-band images consistent with gravitational in- 
teractions. Signs of optical disturbance are also observed in three 
adjacent state A spirals, CGCG 097-092-1, CGCG 097- 138 and 
[BO85]074 showing that tidal interactions, if important, could be 
affecting spirals in more than one evolutionary state. The high rate 
of disturbance in the optical morphology in this region suggests the 
spirals are subject to strong interactions possibly related to their lo- 
cation at the intersection of the NW sub-cluster and the filament 
oriented toward the Coma cluster. 

The state D spirals all have a red colour, except CGCG 097- 
093 which is exceptionally blue (g-i — 0.42) and has a disturbed 
stellar disk. The observed red colours and high Hi deficiency are 
both expected results of a spiral making a transit through a high 
density ICM core. An example of this is found in Virgo, i.e. the 
galaxy NGC4568 (Boselli & Gavazzi 2006a). The calculations for 
the NW sub-cluster (section |5~2| l suggests that a core transit would 
be required to cause the Hi deficiencies observed in the state D 
spirals under the ram pressure stripping scenario, although the 
ICM core transited may have been the SE subcluster. The tendency 
for the state D spirals further from the cluster centre to be redder 
and have higher velocities may be indicating that these spirals are 
a population of backsplash galaxies (Ellingson 20041, which were 
ram pressure stripped as they passed through the SE sub-cluster 
ICM core with a significant velocity component directed to the 
west. 



6 CONCLUDING REMARKS 

We have presented medium-resolution VLA data for a field located 
to the NW of A 1367's centre, and some higher-resolution, but less 
sensitive, VLA data for a field containing BIG. By using a cata- 
logue of the brightest late-type members of the cluster, combining 
AGES and SDSS data, we have identified four evolutionary states 
with shared characteristics in terms of their combination of H I con- 
tent and colour. 

The overall picture from the HI deficiencies as a function of 
radial distance and calculations of ram pressure for the NW clus- 
ter (section [5T2^ are both consistent with moderate levels of ram 
pressure stripping. However the presence of spirals with a range of 



evolutionary states in projected and velocity space, particularly the 
spirals with more advanced evolutionary states (state C and D spi- 
rals), suggests an evolutionary process or processes operating on 
relatively local scales. 

Consideration of the membership of the evolutionary states 
and individual spirals presents mixed, and some times contradic- 
tory, indicators of varying reliability for the dominant ISM removal 
process(es). The lack of correlation between Hi deficiency and 
proximity to the NW sub-cluster core for the spirals likely to be af- 
fected by it, the presence of starbusts, the presence of galaxies with 
H I offset directions inconsistent with ram pressure, together with 
the elevated number of morphologically disturbed red and moder- 
ately Hi deficient spirals (state C) possibly linked to interactions 
during infall, all suggest tidal effects. On the other hand the pres- 
ence of the highly H I deficient (state D) spirals and the possible 
case of gas falling back into the galaxy's potential following an 
ISM core transit (CGCG 097-072) support the ram pressure strip- 
ping model. 

While the analysis above hints that different ISM removal 
mechanisms may be active in this cluster further observations are 
required to discriminate between the hydrodynamic and tidal mech- 
anisms. If anything, this paper shows that the picture in Abell 1367 
is somewhat confusing, probably the result of it being less evolved 
than Coma and Virgo. 
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Appendices 

APPENDIX A: BIG 

The Blue Infalling Group is an interesting case of a compact group 
that has been suggested to be infalling toward the core o f the SE 
sub-cluster with a relative velocity of ~2000 km s _1 iCortese 



|et al.|20 06). As well as the two principal galaxies CGCG 097-1 14 
and CGCG 097-125, BIG contains an extraordinarily high concen- 
tration of dwarf galaxies with high star formation rates and extra- 
galactic H II regions ( Sakai et al. 2002 , Gavazzi et al. 2003a). Deep 
Ha imaging has revealed extensive gas streamers interpreted as 
resulting from a strong tidal interaction between members of the 
group (Cortese et al. 2006 ). H I in BIG has previously been imaged 
using the Westerbork Synthesis Radio Telescope (WSRT) by |Sakai| 
|etaLl ( [2002) . 

Our observations have higher spatial resolution that clearly 
show the offset between the highest density Hi and the optical 
galaxy (Figure [12}. Weak Hi emission was detected by us at the 
position of a knot 15 arcsec (6 kpc) east of knot K2a, with peak 
Hi emission at 8190 km s _1 or ~70 km s _1 greater than K2a's 
optical velocity (Figure \T2\. However, we did not pick up the large 
continuous region of low density H I extending west of CGCG 097- 
125 and running from below K2a north to DW2 seen with the 
WSRT ( Sakai et al. 2002 1, which is probably due to our lower sen- 
sitivity. We were also unable to confirm in our robust cube the 
H I to the west of the optical position of the irregular blue galaxy 
CGCG 097-1 14, observed with the WSRT, although weak emission 
was found after applying natural weighting. The weak detection of 
CGCG 097- 114 implies a large Hi deficiency which is confirmed 



by the Arecibo determined H I mass of 3 x 10 M© l Cortese et al. 
|2006l >. 

Regarding the relative location of BIG with respect to the clus- 
ter centre, this is still matter of debate (Sakai et al. 2002; Gavazzi et 
al. 2003a; Cortese et al. 2006) . The calculation discussed in section 
5.2 indicates that, if the group is within ~ 0.5 Mpc from the cluster 
centre, we should expect evidence of strong ram pressure stripping. 
This, however, would seem to be contradicted by the presence of 
substantial amounts of extended H I in the velocity range from 8000 
to 8500 km s _1 , suggesting either that BIG lies in the background 
of A1367 or it has just started its infall into the cluster center. 
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5500 6000 

Wavelength (1) 



Figure Bl. A 1367 [GP82]1227 emission line spectrum: SPM emission line 
spectrum showing Hp, [OIII] (4958.9), [OIII] (5006.8) [NII](A 6584) and 
H a . 



APPENDIX C: CHANNEL MAPS OF THE VLA 
DETECTIONS 

The channel maps for all VLA detections are available as supple- 
mentary online material. Below as an example of the online mate- 
rial are the channel maps for CGCG 097-079. 



APPENDIX B: ABELL 1367[GP82]1227 



We obtained a 
A1367[GP82]1227 
7225 A, in May 2007 (Figure 



low resolution (4.8 A) spectrum of 
coveri ng the wavelength range 4630- 



Bl 



The observation was carried out 



with the 2.12-m telescope at the San Pedro Martir Observatory 
(SPM) in Mexico, using a Boiler & Chivens spectrograph, a 600 
1/mm grating and a 2K Thomson CCD. This configuration gave a 
dispersion of 1.3 A per pixel. The slit was placed in an East-West 
direction passing through the central region of the galaxy. Four 
20 minute exposures were taken and the spectra were averaged 
after cosmic ray cleaning, extraction and wavelength calibration. A 
wavelength solution using about 30 He+Ar+Ne lamp lines, gave an 
uncertainly of ~0.18 A rms. Sky transparency was good, with the 
seeing slightly above 1 arcsec. Due to the faintness of the galaxy, 
the signal to noise obtained was only about 5 in the range 5500 - 
6500 A. 

The spectrum of A 1367 [GP82] 1227 shows the emission lines 
of H/j (A rcst = 4861.3), [OIII] (4958.9), [OIII] (5006.8) and H Q 
(6562.8) at 4963.8, 5063.2, 51 12.3 and 6700.4 A, respectively. This 
gives an emission-line radial velocity of 6307 km s _1 , which is 
in good agreement with our Hi velocity of 6239 km s~ (see 
Table |2j. These emission lines indicate active star formation in 
A 1367 [GP82]1227. The [OIII] (A 5007)/^ and [Nil] (A 6584)/H Q 
ratios are 2.64 and 0.03, respectively (note that [Nil] is barely vis- 
ible in the spectrum) putting this galaxy completely in the locus of 
star forming galaxies in a BPT diagram (Bald win et al.|1981) . It 
was not possible to resolve the Ha regions individually with the 
SPM 2.12m. 
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Figure CI. Channel maps for CGCG 097-079. The velocity in km s _1 is shown in the top righthand corner of each frame. Contour levels are at -2, 2, 3, 4, 
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Figure Cl. - continued 
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Figure CI. - continued 



